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The present invention provides an optimal procedure for 
determining in the time domain equalizer coefficients for an 
equalizer, where the equalizer compensates for the effects of 
the communications channel on the transmitted signal. A unit 
pulse is transmitted over the communications channel, and a 
channel impulse response is estimated from the received sig- 
nal. A cost function establishes a mean-square error associ- 
ated with the estimated channel impulse response as compared 
to a desired impulse response signal. The value of the cost 
function varies based upon an offset value between the esti- 
mated channel impulse response and an equalized channel im- 
pulse response. Values of the cost function are determined for 
different offsets, and the offset that produces the smallest cost 
function value (corresponding to the minimum mean-square 
error) is selected. The optimal equalizer coefficients are then 
calculated using the selected offset and the established cost 
function. In a preferred embodiment, the invention is applied 
to multicarrier modulation systems. Moreover, the cost func- 
tion is defined as a krakovian function that depends on the 
offset value. That krakovian function includes a first krako- 
vian which uses the channel impulse response and a second 
krakovian which uses the autocorrelation of the channel im- 
pulse response. Using the first and second krakovians and the 
selected offset corresponding to the minimal value of the es- 
tablished cost funciton, the optimal equalizer coefficients are 
determined in a straightforward, one step calculation. 
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AN EQUALIZER FOR USE IN MULTICARRIER 
MODULATION SYSTEMS 

FIELD OF THE INVENTION 

The present invention relates to equalization, arid more particularly, for 
equalizers used in multicarrier modulation systems. 

BACKGROUND AND SUMMARY OF THE INVENTION 

Asymmetric high speed Digital Subscriber Line (ADSL) and Very high 
speed Digital Subscriber Line (VDSL) are examples of modern communication systems 
which permit transmission of data over communication lines at very high rates, e.g., up 
to 52 Mbits/s. One example of a high speed data application is video on demand. The 
transmission of high-speed data over band-limited channels may be accomplished 
through the use of DMT-based, digital communications systems. DMT modems follow 
the general principle of transmitting data by dividing it into several interleaved bit 
streams, and using these bit streams to modulate several carriers. 

A significant limitation in high data rate communication systems is 
intersymbol interference (ISI). One way to compensate for ISI in a DMT system is to 
add a cyclic prefix (guard time) to the beginning of each transmitted DMT symbol. 
Unfortunately, while increasing the length of prefixes reduces ISI, it also decreases the 
effective data rate. Another approach to combating ISI is to employ an equalizer at the 
receiver, but many equalizers require considerable and ongoing computation 
"overhead." There is a need, therefore, to keep the cyclic prefix length as small as 
possible and still compensate for ISI using an equalization technique that does not 
require complicated and ongoing computations. 
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For purposes of this application, equalization is the process of correcting 
or compensating for the ISI caused by the communications channel. In practical 
communications, the frequency response of the channel is not known. Accordingly, an 
equalizer is designed using numerous parameters that must be adjusted on the basis of 
measurements of a channel's signal-affecting characteristics. 

A transversal filter having a delay line with taps spaced by T-seconds 
(where T is the sampling interval and/, = \/T is the sampling rate at the receiver) is a 
common choice for an equalizer. The outputs of the filter taps are multiplied by a filter 
coefficient, summed, and fed to a decision device. The tap coefficients correspond to 
the channel parameters and represent to some extent a "model" of the channel. If the 
coefficients are properly selected, the equalizer significantly attenuates or practically 
removes the interference from symbols that are adjacent in time to the desired symbol. 
The selection of the coefficient values is typically based on minimization of either peak 
distortion or mean-squared distortion. The coefficients are selected so that the equalizer 
output is forced towards zero at //sample points on either side of the desired pulse. 

There are two general types of automatic equalization. The first 
equalization approach transmits a training sequence that is compared at the receiver 
with a locally-generated or otherwise known training sequence. The differences 
between the two sequences are used to set the equalizer filter coefficients. In the second 
approach, commonly referred to as adaptive equalization, the coefficients are 
continually and automatically adjusted directly from the transmitted data. A drawback 
of adaptive equalization is its computational "cost" involved in continually updating the 
filter coefficients so that the channel model is adapted to the current conditions of the 
channel. Equalizer coefficient computational methods employing adaptive algorithms 
such as the least mean square (LMS) or recursive least square (RLS) are 
computationally expensive. At a very high sampling rate for high data rate 
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communications system, this kind of continual LMS updating of the equalizer filter 
coefficients is particularly expensive. 

The present invention seeks to reduce both computational cost and cyclic 
extension, but at the same time effectively equalize received signals to compensate for 
ISI using a short length equalizer. For many communication applications, such as 
VDSL and ADSL, the communications channel does not change that much during a 
particular transmission. In this context of relative "constant" channel conditions, it is 
satisfactory to estimate the channel coefficients during an initial training sequence 
before the message is transmitted without thereafter adjusting the filter coefficient 
values during the transmission. Of course, it may be necessary to monitor channel 
quality during the transmission, and if channel quality decreases below a certain 
threshold, the equalizer coefficients could then be updated. 

The present invention provides an optimal procedure for determining in 
the time domain equalizer coefficients for an equalizer, where the equalizer 
compensates received signals that are distorted by passing through the channel . A unit 
pulse is transmitted over the communications channel, and a channel impulse response 
is estimated from the received signal. A cost function establishes a mean-square error 
associated with the unequalized channel impulse response as compared to a desired 
impulse response signal. The value of the cost function varies based upon an offset 
value between the estimated channel impulse response and an equalized channel 
impulse response. Values of the cost function are determined for different offsets, and 
the offset that produces the smallest cost function value (corresponding to the minimum 
mean-square error) is selected. The optimal equalizer coefficients are then calculated 
using the selected offset and the established cost function. 

In a preferred embodiment, the invention is applied to multicarrier 
modulation systems. Moreover, the cost function is defined as a krakovian function that 
depends on the offset value. That krakovian function includes a first krakovian which 
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uses the channel impulse response and a second krakovian which uses the 
autocorrelation of the channel impulse response. Using the first and second krakovians 
and the selected offset corresponding to the minimal value of the established cost 
function, the optimal equalizer coefficients are determined in a straightforward, one step 
calculation. 

A better understanding of the features and advantages of the present 
invention will be obtained by reference to the following detailed description of the 
invention and the accompanying drawings which set forth an example embodiment in 
which the principles of the invention are utilized. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram illustrating a basic prior art multicarrier 
transmitter, communication channel, and receiver; 

Fig. 2 is a diagram conceptually illustrating frequency domain discrete 
multi-tone modulation (DMT); 

Fig. 3 is a block diagram depicting the principle components of a DMT 
communications system in which the present invention may be advantageously 
employed; 

Fig. 4 is a more detailed block diagram illustrating a time domain 
equalizer in accordance with the present invention; 

Fig. 5 is a flowchart diagram illustrating procedures for setting the time 
domain equalizer coefficients in accordance with an example embodiment of the 
present invention; 
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Fig. 6 is a graph showing the normalized channel impulse response as 
well as the same channel impulse response after equalization in accordance with the 
present invention as a function of time/sample number; 

Fig. 7 is a graph showing an enlarged portion of Fig. 6; 

Fig. 8 is a graph illustrating the cost function (J) in accordance with the 
present invention for various offset values (hq); 

Fig. 9 is a graph showing the spectra of the communications channel and 
equalizer in accordance with the present invention, and the resulting spectrum; and 

Fig. 10 is a graph showing signal-to-noise (SNR) ratios of an equalizer 
used in a DMT receiver in accordance with the present invention for various DMT 
cyclic extension values. 

DETAILED DESCRIPTION OF THE INVENTION 

In the following description, for purposes of explanation and not 
limitation, specific details are set forth, such as particular modulation/demodulation 
techniques, applications, protocols, etc. in order to provide a thorough understanding of 
the present invention. For example, the present invention is described in the context of 
a multicarrier modulation digital communications system in which krakovian algebra is 
used to facilitate calculation of optimal equalizer coefficients. However, it will be 
apparent to one skilled in the art that the present invention may be practiced in other 
embodiments that depart from such specific details. In other instances, detailed 
descriptions of well-known methods, protocols, devices, and circuits are omitted so as 
not to obscure the description of the present invention with unnecessary detail. 
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Multicarrier Modulation System 

Fig. 1 shows an example Discrete Multitone (DMT) communication 
system in which the present invention may be employed. In a Discrete Multitone 
system, in general, the input bit stream is first serial-to-parallel converted. The parallel 
output is then grouped into N groups of bits corresponding to the number of bits per 
symbol. Portions of bits are allocated to each DMT carrier. 

More specifically, transmitter 10 includes a serial-to-parallel converter 14, 
a multicarrier modulator 16, and a pretransmit processor 18. Receiver 12 includes a 
post channel processor 20, a multicarrier demodulator 22, and a parallel-to-serial 
converter 24. The transmitter and receiver are linked in this example by a digital 
subscriber line (DSL) or other form of communication channel 26. Serial input data at a 
rate of b total I T symb bits per second are grouped by converter 14 into blocks of b, ota , 
bits for each multicarrier symbol, with a symbol period of T symh . The b tolal bits in each 
multicarrier symbol are used to modulate N separate carriers in modulator 1 6 with b t 
bits modulating the i-th carrier. 

Preferably, an Inverse Discrete Fourier Transform (IDFT) is used as a 
DMT modulator to generate N s time-domain samples of a transmit signal for each 
block of b total bits, where N s is preferably equal to twice the number of modulated 
carriers, i.e., 2N. The corresponding multicarrier demodulator performs a Discrete 
Fourier Transform (DFT), where bj bits are recovered from the sub-carrier. As 
depicted in Fig. 2, the carriers or subchannels in a DMT system are spaced \/T Hz 
apart across N/T Hz of the frequency band. More detailed discussion of the principles 
of multicarrier transmission and reception is given by J. A. C. Bingham in "Multicarrier 
Modulation for Data Transmission: An Idea Whose Time Has Come", IEEE 
Communications Magazine, Volume 28, Number 5, pp. 5-14, May 1990. 
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In a Discrete Multitone Modulation system, {jf 0 X l9 .. m , X N _ } } are the 
original, complex, input data symbols, {s n } is the modulated data sequence (before 
cyclic prefix), {/?„ } is the discrete-time channel impulse response (CIR), [rj n } is an 
additive noise sequence, \x n } is the received sequence. Adding the cyclic prefix is a 
discrete-time technique used to eliminate ISI in the DMT system. However, as 
described in the background, it is desirable to minimize the length of the prefix in order 
to maximize data throughput assuming the receiver equalizer can effectively remove or 
otherwise compensate for ISL The independent modulating and demodulating vectors 
are the IDFT and the DFT given by the following general relationships: 

2N _ X z&L 

X k= S x n e N * = 0,l,...,2tf-l (1) 
n = 0 



2N-\ 



X 



k 



1 2yv " 1 N 



2N ~ 



-27/-1 (2) 



A DMT system with N channels for transmitting complex data requires a 
DFT size of 2N because of forced Hermite symmetry in the frequency domain. The 
Hermite symmetry is necessary in order to obtain real-valued samples in the time 
domain. In a preferred embodiment, IDFT and DFT are implemented using EFFT and 
FFT algorithms using, for example, an approach disclosed in commonly-assigned 
"Method and Apparatus for Efficient Computation of Discrete Fourier Transform and 
Inverse Discrete Fourier Transform," U.S. patent application Serial No. 08/887,467, 
filed on July 2, 1997. 

Fig. 3 is a block diagram showing the principal operative components of a 
basic DMT transmitter 30 and a basic DMT receiver 32 connected through a channel 
34. Serial input data are grouped into blocks, converted to a parallel form, and approp- 
riately encoded by an encoder 36. Each one of the N active subchannels contains 
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a number of bits allocated to that subchannel. The total number of bits in a DMT 
symbol is therefore: 

N 

B total = lL b k (3) 
k = \ 

Each of the N bit groups is mapped by the symbol into a two-dimensional 
signal constellation. The relation between the number of bits in each block of bits bfc 
and the necessary number of signal points in a constellation Mis given by: 

M k = 2 b " (4) 

The signal points in a constellation are typically arranged in a rectangular 
pattern for even values of and a cross for odd values of b^. Provided that the same 
signal energy is used for all constellations, the distance between neighboring signal 
points is reduced with increasing constellation size resulting in increased BER with the 
same signal-to-noise ratio (SNR). 

The output from the encoder is N complex numbers, one for each group of 
bits, which are then fed to a device that calculates the Inverse Discrete Fourier 
Transform (IDFT). The output is a real sequence that can be considered the 

superimposing of N modulated orthogonal carriers spaced A/apart, Af = — . 

The parallel IFFT output is converted back to a serial data stream by a 
converter 40. The digital modulated data stream is cyclically prefixed, converted to 
analog form by a digital-to-analog converter (DAC) 42, low-pass filtered at 44, and 
passed through a splitter filter and D.C. isolating transformer 46 during pre-transmit 
processing to produce an analog signal transmitted over a transmission channel 34. 

At the receiver, the received analog signal is passed through a D.C. 
isolating transformer and low-pass filter 48, converted to digital form by an analog-to- 
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digital converter (ADC) 50, time domain pre-equalized by a finite impulse response 
(FIR) filter 52 to limit the effective memory of the channel, and stripped of the cyclic 
prefix during post-receive processing in converter 54. The resulting digital signals are 
demodulated by 2N FFT operation 56 and converted to parallel frequency domain 
5 signals. Since the amplitude vs. frequency and the delay vs. frequency responses of the 
channel are not necessarily constant across the entire frequency band, the received 
signal differs from the transmitted signal, and the parallel inputs to the frequency 
domain equalizer (FEQ) 58 differ from the parallel outputs from the encoder 36. It is 
the task of the FEQ 58 to compensate for these differences. The equalized signals can 
)0 then be properly decoded and converted back to a serial form by the decoder 59. 
Ideally, the detected output serial data from the decoder 59 is the same as the input 
serial data to the encoder 36. 

The preferred example embodiment of the present invention employs a 
Time Domain Equalizer (TEQ) 52 which shortens the overall channel impulse response 

15 and therefore reduces ISI. The present invention finds particular advantageous 

application in the training of the time domain equalizer 52. During a training mode of 
operation, a switch 51 is actuated (for example by a general data processor controlling 
the receiver) to connect the output of the analog-to-digital converter 50 to an equalizer 
filter coefficients processor 53 which calculates optimal TEQ coefficient values for a 

20 particular communications channel. Those optimal values are used to set the 

coefficients in equalizer 52. After training, the switch 51 is moved to directly connect 
the sampled output from the ADC to the equalizer 52 for operation in the normal 
equalization mode. 

Fig. 4 illustrates the equalizer 52 in additional detail. In the training 
25 mode, the received sampled signal x n is routed to the equalizer filter coefficients 
processor 53 for calculation of the filter coefficients 

TEQ 0 , TEQ, ^..TEQ^.j , and TEQ A ,_, in accordance with calculation procedures 
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described in detail below. The equalizer 52 is essentially a transversal filter 

corresponding to a delay line with z _I symbol taps, where z _I corresponds to a 
sampling interval T. The tapped signals are weighted or multiplied by a corresponding 
filter coefficient value TEQ in a corresponding multiplier. The multiplier products are 
summed and output as y n . The filter coefficients, which can be expressed as a vector 
TEQ„, are set to minimize negative effects of interference from symbols that are 
adjacent in time to the desired symbol. In general, the criterion for selecting the TEQ 
coefficients is based on minimizing means-square error caused by such interference. 

As described above, a multicarrier transmitter serial-to-parallel converts 
an input data stream to sub-symbols that are distributed to modulate carriers which are 
orthogonal on the symbol duration (IN + v)T, where T denotes the sampling interval. 
Assuming N sub-symbols distributed over 2N carriers and imposing Hermitian 
symmetry, i.e., X" = ^X" N _ k ) *, the transmitted sequence s n in the time domain is 
given by: 



, oo 2N-1 2 ^ k TiO 
1 x — < < „. 2/v 

^=^Z E*r« (5) 

w=-co *=0 



where X" is the sub-symbol conveyed by sub-carrier k during the m-th symbol, and the 
1 

27V 



factor -5- is introduced to normalize the overall transmission power. 



The overall channel is considered to include the physical channel plus the 
filtering performed in the transmit and receive circuitry at both ends of the channel. 
The channel impulse response (CIR) is assumed to be linear. Denoting CIR by h(t), 

and providing cyclic extension consisting of v samples, the received analog signal is 
given by: 
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*(') = ^7 t Y Y X " h {t - nT - m(2N + v)^™ + rfc) 

Z/V m=-oo k=0 n=-v 

= 'S^rE' h(t - nT - m(2N + v)f)e^ + rj(t) (6) 

m=-oo n=— v )fc=0 

where is the additive white gaussian noise (AWGN) with one-sided spectral 
density /V 0 . 

At the receiver, the cyclic prefix samples are discarded, and FFT 
(demodulation) is performed on the IN real-valued samples. Inspecting equation (6), 
the innermost summation corresponds to the Fourier coefficient of the CIR associated 
with the sub-carrier k. Hence, after demodulation, each symbol only needs to be 
corrected by a complex coefficient associated with the channel transfer function that 
includes a phase shift and attenuation. Alternatively, this correction in the time domain 
may be seen as convolution of the demodulated sequence with the composite channel 
impulse response. 

As described earlier, the channel causes ISI which may severely 
deteriorate system performance (e.g., poor symbol decisions). The dispersiveness of the 
channel may be compensated for by adding guard time. For DMT systems, each 
symbol is cyclically extended (the v last samples are repeated at the beginning of the 
block data) so that input sequence "looks" periodic. Again, the disadvantage of cyclic 
extension is reduced throughput efficiency which is proportional to the extension 

2N 

length, i.e., by a factor ^ ^ . For severe-ISI channels, throughput efficiency 

becomes intolerably low. Increasing Nor v increases computational complexity, 
system delay, and memory requirements of the transceiver, and hence the cost of 
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transceiver increases quite significantly. This problem is solved in the present invention 



Mathematical Principles Used to Formulate the Equalizer 

The complete suppression of ISI requires that the equalized signal, y n , be 
5 identical to the transmitted sequence, s n , up to an offset n 0 (defined greater detail 
below): 



where y represents an arbitrary scaling factor. If the original sequence X k to be 
transmit is complex, it is reformatted as a Hermite symmetric sequence. To simplify the 
10 notation, the index "m" which represents the symbol number is dropped. Thereafter, the 
IDFT of the Hermite symmetric sequence is calculated, (which corresponds to 
modulation), producing a real signal which is transmit over the channel. The received 
signal is then demodulated using the DFT. 

The equalizer in accordance with the invention is therefore designed to 
15 determine a vector TEQ„ such that when TEQ„ operates on (or filters) the sampled 
received signal, x n , the result, y n , is (to a great extent) a replica of the originally 
transmitted sequence s„, i.e., 



using a short-length, computationally "inexpensive" equalizer. 



(J) 



y n = TEQ n ®x { 



n 



(8) 



Substituting x„ using equation (6) results in: 




n 



1 2iV-l 

Jfi L lLTEQ n ®h(t- n T-rr(2N+v)T)Y,X\ 



Irgk 



2N 



20 



y n = 



+ In (9) 
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where rj' is filtered external noise, and the symbol ® depicts convolution. 

In order for the equalizer to perfectly suppress ISI, it must satisfy the 

following: 

TEQ n ®h n = r -S{ n - no ) + ^' (lb) 

where 7^' represents the deviation of the overall channel impulse response (CIR) from 
the desirable (ideal) unit pulse S(n) , i.e., the error due to imperfect equalization. The 
goal is to compensate the CIR (h n ) so that ^ is reduced to zero. Therefore, in the 
absence of noise and using an infinite length TEQ n , it is theoretically possible to 
achieve complete suppression of ISI. But in practical applications, the equalizer is of 
finite length. Assuming an offset n' Q associated with unknown delay introduced by the 
channel n 0 , channel model deviation (or channel misfit) 77,;' , and external noise sources 
77,; , the mean error may be minimized towards zero using the Least Squares principle. 

Examining the limited number of equalizer coefficients TEQ n necessary 
to minimize the error ^ , assuming an originally transmit dirac pulse 5(n), we insert 
equation (9) into equation (10) and obtain: 

I cd 2N-\ 2N-\ " 

y » = yu?< Z(r^(«-« 0 )+^')S^ 2 ^ + /7; 

I °° 2/V-l 2 . k n 

= r E Yv»lLXLe 1Nj r< (11) 

m=-oo n=- v k=0 

where Tj' n " depicts the noise sequence due to the imperfectly equalized channel. 
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The filtering operation, as described by equation (9), involves the 
convolution of an infinitely long sequence with a finite one. Using krakovian algebraic 
notation, equation (10) may be rewritten as the following: 



' TEQ 0 ^ 




K\ 


•■ K • 


h 1 








£ m\ 


< 










> — y ■ < 


1 


> + < 








h m\-(K-\) ■ 




" ^m2-(AT-l)j 




0 







where: indices ml -n 0 and ml — > oo, h 0 represents a first maximum amplitude of 
the sampled CIR (before equalization), y - r {0, . . . , 0 , 1, 0 , . . . , 0} represents the desired 
CIR, and e k is error sequence comprising both filtered external noise and noise due to 
imperfect channel equalization referred to hereafter as an observation error s. 

The following is a brief and basic description of krakovian algebra. A 
krakovian is defined as a set of elements (numbers, symbols, or other krakovians) 
arranged in rectangular array of m columns and n rows so that the krakovian has a 
dimension m x n . The elements of a krakovian are indexed according to their position, 
namely by column and row number, krakovian algebra is not the same as matrix 
algebra. In fact, some of the different properties of krakovians, not found in regular 
matrix algebra, are employed in the preferred embodiment to achieve a'significant 
reduction in data processing overhead. 

In krakovian algebra, the product of two krakovian factors a and b when it 
consists of accumulated cross-products of the columns of krakovians a and 6. The 
product clement is written in the product column having the same column number as the 
first factor a and in the product row having the same number as the column of the 
second factor b: 
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Py=Y, a tk b jk (13) 

Thus, krakovian multiplication requires the same number of rows. If the product a • b 
exists, then the product b • a always exists and is equal to the transpose of a • b . This 
definition leads to a laws of commutation and association that are different from the 
commonly known matrix laws. 

A special krakovian, referred to as the diagonal krakovian, includes all 
"1 V on the NW-SE diagonal and "0's" in every other location and is depicted by the 
Greek letter r . When the krakovian r occurs as a second multiplicand, it plays a 
similar role in the krakovian algebra as does unity in the arithmetic of ordinary 
numbers. However, a krakovian r as the first multiplicand causes a krakovian 
transpose. These statements may be expressed as follows: 

a-r = a (14) 

r-a = ra (15) 

where the transpose of a krakovian means as it does in regular matrix algebra: 

i Ta ) iJ =a ji (1.6) 

Changes in the order of multiplying krakovians follow certain rules. The 
krakovian association law states: 

a b o d e- f = a -(^e- rd ■ tc ■ rb)- f (17) 
The krakovian law of dissociation is 

a(bcde)f = aTerd-vcbf (18) 
Krakovians do not fulfill certain matrix algebra relationships, e.g., 
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(ab)c*a-(b-c) (19) 

With this very brief explanation of krakovian algebra, we return to the 
derivation of the optimal equalizer filter coefficients. Equation (12) may be rewritten i 
more compact form in krakovian notation 

TEQ n xa = / + e (20) 

where TEQ n is a vector containing coefficients TEQ 0 ...TEQ K _ X of me time domain 
equalizer 52, a is a krakovian containing samples of the received channel impulse 
response (CIR), / is a krakovian vector showing a "best" sample location for the 
received unit pulse, £ is a vector representing the observation error, and r denotes 
transpose. The degree of the TEQ„ satisfies K < n 0 + m2 since equation (12) 
represents an over-determined system Hereafter, a description of krakovian 
dimensions is omitted when clear from the context. 

The krakovian vector / on the right side of the equation (20) is 
purposefully configured to have only one non-zero element at the offset position n' Q , 
where -n 0 <n' 0 < ml . Vector / will be referred to as the "coordinate vector" because 
of a physical interpretation as the delay introduced by a communication channel. The 
offset n' Q is the sampling time difference between the first maximum peak of the CIR 
and the peak of the equalized CIR. As is shown below, the choice of the offset n' 0 
affects minimum mean-square-error. Hence, n' 0 must be selected to minimize the mean 
square error. 

The mean square error may be viewed as all of the output energy except 
for the energy transmitted in position n' Q , and is established as a cost function /(«o) 
corresponding to the square of the observation error. 
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4"°) = {(n®. ® *.) - y ■ s(» - »;)}* (21) 

The cost 7 is calculated at multiple offsets n' Q . The optimal coefficient vector TEQ„ , 
(i.e., the optimal filter coefficients), may be obtained when the cost function j(n' 0 ) 
reaches a rninimum. If the solution of equation (21) exists in a mean-square sense, the 
krakovian a containing the sampled CIR, h k , must fulfill the following constraint. 

The square of a must be "positive definite" which means that (a) 2 must be capable of 
being inverted. 

The following assumptions and manipulations, using krakovian algebra, 
are useful in efficiently solving equation (21) for minimum mean-square error. 
Equation (20) may be expressed using krakovian notation as: 

TEQ' -Ta' = s (22) 

where TEQ' = r{TEQ„ «}, u = -1, and a' = {a /„ 6 }. This assumption permits 

squaring of both sides since equation (22) provides the error in a "pure" form which 
permits easy calculation of the cost function J 

./ = (TEQ'.ra') 2 = * 2 (23 ) 

Equation (23) may be manipulated using krakovian algebra, resulting in the following: 

s 2 = (TEQ' • ra') 2 = TEQ' -(a') 2 • TEQ' = (TEQ' • rr") 1 (24) 

Given that the krakovian a corresponds to the CIR, the krakovian (of corresponds to 
the autocorrelation (AC) of the CIR. As is well known, the autocorrelation function 
and power spectral density are Fourier transform pairs. The receiver demodulator in a 
DMT system can be used to generate the autocorrelation of the CIR. 
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Taking the square root of equation (24) results in 

TEQ' - tr" = t{w 0 ... Wk _ t w K } (25) 
Recall that the goal is to minimize the square error which can now be expressed as 

K 

£ 2 = 2! w k ~ minimum (26) 

k=0 

5 The Krakovian root r" has the structure given by the following notation: 

r "-{o 3 (2? > 

where r is the upper triangle with positive diagonal entries, 0 depicts zeros below the 
main diagonal, and v[p w) is the last column of r" . Using the above notation, we 
obtain from equation (25) 

10 TEQ n -zr + up=T{w a ...w fC _ ] } (28) 

w=w K (29) 

Equation (28) always has a solution irrespective of the choice of w 0 ,..., w K _ } . 

However, the minimum of mean-square-error expressed in equation (26) is achieved 
when 

15 w 0 = ... = m^_, -0 (30) 

In other words, w 2 is the least possible value of e 1 , where e 2 is equal to the cost 
function J(n' 0 ) for a particular coordinate vector l n . Q which can be expressed as: 



min[e 2 ) = w 2 =j( n ' Q ) (31) 
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Using these properties, the sum of squared observational errors is 
calculated without the use of an adaptive algorithm like the gradient search algorithm. 
Instead, a purely algebraic procedure is employed thereby circumventing the difficulties 
associated with the existence of local minima, and achieving the optimal solution in a 
finite number of steps corresponding to the number of choices of n' 0 . 

Equation (24) yields 

{<*'?= {r"Y (32) 
Combining equation (22) and equation (27) with equation (32) results in 

(r) 2 = (af (33) 
P r = l nb a (34) 

* a =('*) 2 -to 2 =l-to 2 05) 

Substituting equations (34) and (30) into equation (28) and using equation (33), we 
obtain a formula for determining the optimum value of the vector TEQ„ : 

TE Q»=-KV4(WT (36) 

The term (/„ 6 • a) corresponds to the « 0 -th column of the krakovian a , and the 

inverted krakovian ((a) 2 ) ' may be physically interpreted as the inverse of a krakovian 
consisting of the autocorrelation sequence of the channel impulse response. 

From equation (36), the corresponding value of the cost function j(n' Q ) 
can then be expressed as: 
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y(^)=l + M TEQ„ •(/„.•«) (37) 

or equivalently 

J{^)=l-u 2 (l nb -a)-((a) 2 y -(l nb -a) (38) 

In this form, the cost function J{n 9 0 ^ may be computed and its minimum determined 
5 without any knowledge of the equalizer filter coefficients TEQ W . Since determining 
the optimum equalizer coefficients TEQ W is the end goal, the traditional "fitting" 
procedure is reversed in the present invention. Instead of searching for the minimum 
mean-square-error by adaptively adjusting filter coefficients, (e.g., using recursive 
LMS), the optimum equalizer coefficients are, determined in a single step as described 
10 below after the search for minimum mean-square-error is accomplished using 
equation (38). 

To determine the minimum mean-square error, the cost function J{n' 0 ) in 
equation (38) is determined for discrete choices of offset n' 0 position. A good initial 
guess about the region of search for the minimum of cost function can be obtained by 
15 the analysis of available communication channels. This cost function determination 
procedure is purely algebraic and discrete, therefore the cost function need not be 
differentiable. In fact, the cost function j(n r 0 ) is non-differentiable in respect to the 
offset n' 0 . Thus, the only necessary assumption is the non-singularity of the inverse 

krakovian ((#) 2 ) ■ The optimum short-length equalizer coefficients TEQ n may be 
20 determined in a single step from equation (36) using the column vector • a obtained 
in equation (38) for which the minimum value of J{n r ^) was obtained in equation (36), 
keeping in mind the krakovian (a 2 ) is a constant. 
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Equalizer Coefficient Setting Procedure 

Having demonstrated the underlying mathematical principles of the 
present invention, a routine (block 60) for determining the optimal equalizer 
coefficients in accordance with an example implementation of the invention is now 
described in conjunction with the flowchart in Fig. 5. Switch 51 is connected to the 
equalizer filter coefficients processor 53 when training the equalizer 52. A unit or dirac 
pulse is transmitted over the channel and received in the receiver (block 62). The 
equalizer filter coefficients processor 53 processes the received samples x n and 
estimates a channel impulse response (CIR) of the communications channel (block 64). 
Using the channel impulse response, the equalizer filter coefficients processor 53 
estimates the autocorrelation function (ACF) associated with the just-estimated channel 

impulse response as a krakovian (a) 2 (block 66). The autocorrelation function 
krakovian is then inverted taking advantage of the Toeplitz structure of the krakovian 

(a) 2 to generate ((«) 2 ) (block 68). Using the estimate of the CIR, the following 
krakovian vector is formed: (/„ & -«) = r(/7„. +| ,h n ^ ,h nb+K] ) (block 69). 

The next task is to calculate the cost function J as a function of the 
channel impulse response offset n' 0 using equation (38) (block 70). Various channel 
impulse response offset values are inserted into block 69 in order to determine that 
offset n' Q for which the cost function J(n' 0 ) is the smallest, i.e., minimum mean-square 
error (block 72). The optimum equalizer filter coefficients TEQ„ are then determined 
using equation (36) (block 74). Once calculated, the equalizer filter coefficients 
processor 53 then sets the equalizer filter coefficients TEQ 0 , TEQ S , ... TEQ K _ 2 , and 
te Qk-\ to the calculated optimal values (block 76). After the coefficients are set, 
switch 5 1 is then moved to the normal position, and the equalizer 52 processes received 
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signal samples x n using the set TEQ filter coefficients. Equalizer output samples y n 
are then compensated for the channel-induced ISI (block 78). 

Performance of an Equalizer Implemented Using the Invention 

Figs. 6-9 are now described to illustrate the performance of an example 
implementation of the present invention by way of simulated tests. A one kilometer 
telephone loop was chosen as an example communications channel. A unit impulse was 
transmitted over that telephone loop. Fig. 6 graphs normalized amplitude on the 
vertical axis versus time (indicated as received sample number) on the horizontal axis. 
The channel impulse response (CIR) received corresponding to the transmitted unit 
pulse is shown as an oscillating waveform (indicated as a dashed line) having multiple 
"pulses" which obviously adversely interferes with adjacent symbols. In contrast, the 
same channel impulse response after equalization (shown as a solid line) includes only 
one very narrow spike shortly after the first maximum pulse of the unequalized channel 
impulse response. Accordingly, the equalized channel impulse response resembles very 
much the originally transmitted unit pulse with minimal overlap onto adjacent symbols. 

Fig. 7 enlarges the portion of the graph in Fig. 6 around the first 
maximum of the CIR from a reference chosen at sample zero. It is easy to see that the 
channel impulse response offset n' 0 , defined as the offset between the first maximum of 
the channel impulse response and the point at which the cost function J(n^ is 
minimized, corresponds to one sampling interval 

Fig. 8 illustrates this point somewhat differently by graphing the cost 
function J(n' Q ) on the vertical axis as a function of the offset n' 0 . The first maximum 
of the unequalized channel response is at time sample zero. By substituting different 
offset n' 0 values, e.g., -3, -2, 1, 2, 3, etc., and calculating the cost function, the graph in 
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Fig. 8 clearly shows that the minimum cost function value is at an offset n r 0 equal to 
one sample time from sample zero. 

Fig. 9 illustrates the overall performance of this example equalizer. The 
spectrum of the one kilometer telephone loop channel is illustrated as the dashed dot 
5 line. The equalizer's spectrum, is shown as the dashed line and substantially mirrors the 
spectrum, of the communications channel. The resulting spectrum equalized in 
accordance with the present invention (shown as a solid line) is for the most part flat 
which corresponds to optimal equalization. 

Comparing the original channel impulse response with the equalized 
10 channel impulse response, it is clear that the present invention offers a substantially 
improved signal-to-noise ratio performance with respect to intersymbol interference. 
Fig. 10 illustrates, for this example in a DMT environment, SNR performance in the 
frequency domain due to the uncancelled, oscillating "tail" of the channel impulse 
response, (which causes intersymbol interference), for various length DMT cyclic 
15 extensions. As is evident from Fig. 10, there is no dramatic difference in SNR even 
though different cyclic extension length values are used (64, 32, and 16). Accordingly, 
using the present invention, shorter cyclic extensions may be used resulting in greater 
throughput efficiency. 

Optimal equalization is achieved in the present invention using a short 
20 length, non-recursive equalizer that has a number of advantages over gradient-based, 
recursive equalizer algorithms. One well known drawback of gradient, recursive 
algorithms is their susceptibility to local minima. The present invention avoids this 
problem using the algebraic principle of minimizing means square error. Instead of 
differentiating the cost function, the cost function's value is explicitly monitored for 
25 discrete choices of the channel impulse response offset vector until a minimum is found. 
Thereafter, the optimum equalizer coefficients of a short length equalizer are readily 
calculated in non-recursive fashion and with relatively small computation cost. 
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While the invention has been described in connection with practical and 
preferred embodiments, it is to be understood that the invention is not to be limited to 
the disclosed embodiments, but on the contrary, is intended to cover various 
modifications and equivalent arrangements included within the spirit and scope of the 
appended claims. 
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WHAT IS CLAIMED IS: 

1 . In a receiver that receives signals transmitted over a communications 
channel where the received signals are compensated in an equalizer for effects of the 
communications channel on the transmitted signal, a method, comprising the steps of: 

estimating a channel impulse response for the communications channel; 
estimating an autocorrelation function associated with the channel impulse 
response; and 

calculating coefficients for the equalizer from the estimated channel impulse 
response and the estimated autocorrelation function. 

2. The method in claim 1, further comprising: 
determining an error associated with imperfect equalization; 
determining a cost function associated with the error; and 
calculating the equalizer coefficients using the determined cost function. 

3. The method in claim 2, wherein the cost function defines a square error 
associated with an equalized version of the channel impulse response compared to a 
desired channel impulse response. 

4. The method in claim 3, wherein a value of the cost function varies 
depending on an offset of the estimated channel impulse response and the equalized 
version of the channel impulse response. 

5. The method in claim 4, further comprising: 
determining the value of the cost function for different offsets, and 
selecting the offset that produces the smallest cost function value, 
wherein the equalizer coefficients are calculated using the selected offset. 



WO 99/50969 



PCT/SE99/00468 



26 

6. The method in claim 1, wherein the channel impulse response and 
autocorrelation are determined using an initial training signal including a unit pulse 
transmitted over the communications to the receiver. 

7. The method in claim 1, wherein the calculating step employs krakovian 
5 algebraic relationships. 

8. A receiver that receives signals transmitted over a communications 
channel comprising: 

an equalizer that compensates the received signal for effects of the 
communications channel on the transmitted signal; and 
10 a processor that determines parameters of the equalizer associated with the 

communications channel based on krakovian algebra. 

9. The receiver in claim 8, wherein the equalizer is a equalizer implemented 
as a finite impulse response (FIR) transversal filter having multiple taps and 
corresponding equalizer coefficients. 

15 10. The receiver in claim 9, wherein the equalizer includes a training mode, 

during which optimal equalizer coefficients are calculated and set by the processor, and 
a normal operation mode, during which the set coefficients are used to weight the 
signals at the corresponding filter tap, the weighted signals being summed to produce an 
equalized signal. 

20 11. The receiver in claim 9, wherein the receiver is a discrete multi-tone 

(DMT) receiver. 

12. The receiver in claim 9, wherein the processor defines a first krakovian 
using a channel impulse response of the communications channel, defines a second 
krakovian using an autocorrelation function associated with the channel impulse 
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response, and calculates the optimal equalizer coefficients using the first and second 
krakovians. 

13. The receiver in claim 12, wherein the processor inverts the second 
krakovian to calculate the optimal equalizer coefficients. 

14. The receiver in claim 9, wherein the processor determines a squared error 
between the estimated channel impulse response and a desired impulse response and 
calculates the equalizer coefficients minimizing the squared error. 

15. The receiver in claim 14, wherein a value of the squared error varies 
depending on an offset of the estimated channel impulse response and an equalized 
version of the channel impulse response. 

16. The receiver in claim 15, wherein the processor determines the value of 
the squared error for different offsets, selects the offset that produces the smallest 
squared error value, and calculates the equalizer coefficients using the selected offset. 

1 7. An equalization method for equalizing received signals transmitted over a 
communications channel for the duration of a transmission using equalization 
parameters, comprising the steps of: 

estimating a channel impulse response of the communications channel using a 
known signal transmitted in the initial portion of the transmission; 

determining an offset value between the estimated channel impulse response and 
an equalized channel impulse response; and 

determining the equalization parameters using the determined offset value. 

18. The method in claim 17, further comprising: 

manipulating the channel impulse response to determine a corresponding 
autocorrelation function; 

assembling the autocorrelation function as a krakovian; 
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inverting the krakovian; 

determining optimal values for the channel model parameters using the inverted 
krakovian. 

19. The method in claim 17, wherein the equalization method is employed in 
a multicarrier transmission system and the manipulating step includes: 

taking the Fourier transform of the received channel impulse response, squaring 
the magnitude of the Fourier transform result, and taking the inverse Fourier transform 
of the squared magnitude to compute the autocorrelation function. 

20. The method in claim 1 7, wherein the channel model parameters include 
equalization coefficients, the method further comprising: 

determining a cost function J according to the following: 

j{(TEQ„<S>h n )-yS(n-n 0 )} 2 
where TEQ n corresponds to the equalization coefficients, h n depicts the channel 

impulse response, <8> corresponds to convolution, y is a scaling factor, and 6(n - n 0 ) 
corresponds to a unit pulse. 

2 1 . The method in claim 20, wherein the cost function J is formulated as a 
krakovian function that depends on the offset value including a first krakovian 
associated with the channel impulse response and a second krakovian associated with 
the autocorrelation function. 

22. The method in claim 2 1 , further comprising: 

solving the krakovian function for different offset values in order to determine 
the minimal cost function, and 

determining the equalization coefficients TEQ„ using the offset value that 

produced the minimal cost function. 
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(57) Abstract 

The present invention provides an optimal procedure for 
determining in the time domain equalizer coefficients for an 
equalizer, where the equalizer compensates for the effects of 
the communications channel on the transmitted signal. A unit 
pulse is transmitted over the communications channel, and a 
channel impulse response is estimated from the received sig- 
nal. A cost function establishes a mean— square error associ- 
ated with the estimated channel impulse response as compared 
to a desired impulse response signal. The value of the cost 
function varies based upon an offset value between the esti- 
mated channel impulse response and an equalized channel im- 
pulse response. Values of the cost function are determined for 
different offsets, and the offset that produces the smallest cost 
function value (corresponding to the minimum mean-square 
error) is selected. The optimal equalizer coefficients are then 
calculated using the selected offset and the established cost 
function. In a preferred embodiment, the invention is applied 
to multicarrier modulation systems. Moreover, the cost func- 
tion is defined as a krakovian function that depends on the 
offset value. That krakovian function includes a first krako- 
vian which uses the channel impulse response and a second 
krakovian which uses the autocorrelation of the channel im- 
pulse response. Using the first and second krakovians and the 
selected offset corresponding to the minimal value of the es- 
tablished cost funciton, the optimal equalizer coefficients are 
determined in a straightforward, one step calculation. 
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